Quantum Dark Soliton: 
Non-Perturbative Diffusion of Phase and Position 
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The dark soliton solution of the Gross-Pitaevskii equation 
in one dimension has two parameters that do not change the 
energy of the solution: the global phase of the condensate 
wave function and the position of the soliton. These de- 
generacies appear in the Bogoliubov theory as Bogoliubov 
modes with zero frequencies and zero norms. These "zero 
modes" cannot be quantized as the usual Bogoliubov quasi- 
particle harmonic oscillators. They must be treated in a non- 
perturbative way. In this paper I develop non-perturbative 
theory of zero modes. This theory provides non-perturbative 
description of quantum phase diffusion and quantum diffusion 
of soliton position. An initially well localized wave packet for 
soliton position is predicted to disperse beyond the width of 
the soliton. 

PACS numbers: 03.75.Fi, 05.30.Jp 



I. INTRODUCTION 

A dark soliton in a quasi one dimensional (ID) atomic 
Bose Einstein condensate (BEC) is a particle-like solution 
of the Gross-Pitaevskii equation [1]. Dark solitons were 
observed in two experiments [2,3]. A classical soliton in 
a quasi- ID condensate in a harmonic trap behaves like a 
(negative mass) particle in a harmonic potential with a 
frequency equal to trap frequency divided by v2 [1] . This 
results in a (negative frequency) anomalous mode in the 
spectrum of the Bogoliubov theory. This mode describes 
small fluctuations of the soliton around the center of the 
trap. A soliton wave packet which is initially localized in 
the center of the trap is going to disperse [7] . The width 
of the wave packet grows until it becomes comparable to 
the width of the soliton - the healing length. In Rcfs. [7] 
this dispersion was estimated to happen, for reasonable 
experimental parameters, on a time scale of 10ms. This 
quantum "instability" , present even at zero temperature, 
adds to the list of more classical decay channels [8,9]. 

The calculations in Ref. [7] were done within Bogoli- 
ubov theory which is a linearized theory of small quan- 
tum fluctuations around a classical soliton localized in 
the center of the trap. This theory breaks down when 
fluctuations grow large because for large fluctuations one 
cannot neglect interactions between Bogoliubov modes. 
The perturbative theory breaks down when the width 
of the soliton wave packet becomes comparable with the 
healing length i.e. after around 10ms from soliton cre- 
ation. To extend the soliton diffusion beyond this point 
one has to treat the diffusion in a non-perturbative way. 



In this paper I develop non-perturbative theory of soli- 
ton diffusion that includes both position and phase fluc- 
tuations. To avoid some technical problems, and also to 
afford more pedagogical presentation, I consider a soli- 
ton in uniform ID condensate confined to a box of finite 
length. The finite length of the box makes quantum and 
thermal depletion from the ID condensate finite. In fact 
finite quasi-lD condensates were observed in experiments 
[4,2] in harmonic traps. A quasi-lD condensate in the 
Thomas-Fermi regime can be considered locally uniform 
on lengthscales comparable to the soliton width. What 
is more atom chips [5] make possible confinement of cold 
atoms in a quasi-lD box. 
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xl and xr are 

placed six healing lengths from the left and right walls of 
the box. 



FIG. 1. The function f (x) = ^?#L. 



The dark soliton [1] is a stationary solution of the 
Gross-Pitaevskii equation in ID 



iUdtcj) = -—dl<\> + g\cb\ 2 (b - n<t>. 
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Here m is atomic mass and g is effective ID interaction 
strength. The system is placed in a box by imposing the 
boundary conditions (f>(x = 0) = and <j){x — I) at the 
walls of the box. The stationary dark soliton is 
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see Figure 1. Here p is (linear) density of the condensate, 
£ = h/^/mgp is the healing length, c = y/ pg/m is the 
speed of sound, and 9 is arbitrary global phase. I made a 
convenient choice of p = gp. I assume that the width of 
the soliton is much less than the size of the box, £ -C I, 
and that the position of the soliton q is at safe distance 
of a few healing lengths from the walls. 

This stationary solution is degenerate with respect to 
the soliton position q and to the global phase 9 of the 
condensate wave function. 



III. SMALL FLUCTUATIONS AROUND DARK 
SOLITON 

The Gross-Pitaevskii equation (1) can be linearized in 
small fluctuations 5<j>(t, x) around the stationary classical 
background (2): 
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Here the linear differential operator is 

r= (+H GP +. g\6 (x) | 2 +g<t>l{x) 

\ -#5 2 (z) -n GP - g \Mx)\ 2 



H GP = -^dl+ g \Mx)f 



The right eigenmodes of the non-hermitian C are solu- 
tions of the Bogoliubov-de Gennes equations 
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Every right eigenmode has a corresponding left eigen- 
mode (u*,—v*). The right eigenmodes can be classified 
as phonon modes (e > 0) and zero modes (e = 0). 



The soliton does not scatter phonons but shifts their 
phase: as u k (x) or v k {x) is passing from left to right the 

function + ztanh^^^ is changing phase by Aip k — 
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(4)- 



B. Phonons in a box 



In a box the wavevector k is quantized. In figure 1 we 
can see a half-antikink at x = 0, a kink at x = q, and 
another half-antikink at x = In this background total 
phaseshift of a phonon passing from the left wall to the 
right wall of the box is — ^A(p k + Af k — \A(p k = and 
the quantization condition is the same as for a quantum 
particle in a box: k n = with n = 1, 2, 3, . . .. Discrete 
Bogoliubov modes are 



u+k n (—x,0) — U-k n (—X,0) , XL < x 
U„(x) = { u +kn (x,q) - u^ kn (x,q) , x L < x < x R 

U +kn (-X, -I) - U-hJ —x,—l) , x R < x 

plus a formula for v n (x) obtained by replacing u with 
v. With proper normalisation phonon modes satisfy the 
orthogonality relation: (u m \u n ) - (v m \v n ) = 5 mn . 



C. Zero modes 

In addition to phonons there are two zero modes with 
e = 0. One originates [6] from the global i7(l) gauge 
invariance 4>e~ %6 — > (j)e^^ e+f ^ broken by the classical so- 
lution (2): 
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A. Phonons 



and the other from the translational invariance q 
broken by the solution (2): 



q + e 



For later convenience I write first continuous spectrum 
e k and cigenfunctions for phonons on the background 
<j>o(x) — e^y/ptanh extending to x = ±oo, 
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v k (x,q) = 1L e kx e ie x 
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The zero modes are orthogonal: (ug\u q ) — (vg\v q ) = 0. 

Unlike phonons both zero frequency modes also have 
zero norms: (u\u) — (v\v) =0. As a result, phonon modes 
together with zero modes do not span the whole Hilbert 
space in the functional space (S<j>,6<j)*). For example, to 
find a coordinate of (5<j>, 8<f>*) in the direction of the zero 
mode (ug,vg) one has to project (Sep, 6(f)*) on an adjoint 
vector (uf,vf), 
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which has unit overlap with the zero mode, (ug d \ue) — 
(vg d \vg) = 1, but is orthogonal to all other modes. 
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In a similar way, the translational mode (u q ,v q ) re- 
quires an adjoint mode (u q ,v q ) normalized so that 
{iiq d \u q ) — (v q d \v q ) = 1, but orthogonal to all other modes 



including (ug d 



In summary, two adjoint modes 



are missing in order to span the whole Hilbert space of 



D. Adjoint modes 

In the case of the dark soliton (2) a vector (u ad ,v ad ) 
adjoint to a zero mode (u, v) turns out to be a solution 
of the inhomogeneous equation [6] 
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with a constant M chosen so that the overlap (u ad \u) — 
(v ad \v) — 1. Eq.(8) warrants that the adjoint vector 
(u ad ,v ad ) is an eigenstate of C 2 with eigenvalue 0. As 
such it is orthogonal to all phonon modes because phonon 
modes are eigenstates of C? with non-zero eigenvalues e\. 

To find the adjoint vector to the gauge mode (5) it is 
good to start form a stationary Gross-Pitaevskii equation 
solved by <fto in Eq.(2) 
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Taking derivative of this equation and its complex con- 
jugate with respect to p gives 



(10) 



Comparing this equation with (8) and (5) and using 
d P V = g gives 
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The normalization condition (ug \ug) — (v# \ve) = 1 re- 
quires that M§ = Here Nq = (0o |0o) is average 
number of atoms in the condensate mode. With this Mg 
one recovers the general formula [6] 
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To get the adjoint vector to the translational mode (6) 
we verify first that 
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Here the envelope function is 



I(x) 



tanh | 
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x < XL 

X L < X < X R 

Xr<x 



(13) 



(14) 



Comparing this equation with Eqs.(8,6) gives 
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The normalization condition (u q d \u q ) — (v q d \v q ) = 1 re- 
quires M q = - ATip/c < 0. 

The overlap between two adjoint modes should vanish 
but it does not: {u ad \u a e d ) - {vf\v% d ) = iR with real 
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This problem can be easily fixed because 



hgcM q M B ' 

the solution of the inhomogeneous equation (8) is not 
unique - we can always add a zero mode of C to the 
solution. Using this freedom I replace the adjoint gauge 
mode as 
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The new adjoint gauge mode has no overlap with the 
adjoint translational mode. 



E. Decomposition of unity 

Putting together phonons and zero modes with their 
adjoint partners results in a decomposition of the unit 
operator as 



(17) 
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This decomposition is a foundation of the Bogoliubov 
theory 



IV. PERTURB ATIVE THEORY 

The Gross-Pitaevskii equation (1) is a classical version 
of a second quantized theory with a Hamiltonian 



H 
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Here i>(x) is a bosonic field operator. In standard per- 
turbative treatment [6] of zero modes the field operator 
is expanded in small quantum fluctuations Sip around 
the classical solution fa with fixed q = qo and 9 = 0: 
ip = 4> + Sip. The fluctuation operator is expanded as 
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Here the subscript o means a mode with q = q a and 
8 = 0. Fluctuating position q — q and phase 9 are as- 
sumed small. P q and Pg are momenta conjugate to q and 
9: [q, P q ] = ih and [6, Pg] = ih. b n is a phonon anni- 
hilation operator: [6 m ,6t] = S mn . These commutation 
relations plus the decomposition of unity (17) give the 
desired [6x/j(x), 5ip'(y)] = 6(x — y). 

In the Bogoliubov theory the Hamiltonian (18) is ex- 
panded to second order in Sip. The linear term van- 
ishes because the classical <j>o is a solution of the Gross- 
Pitaevskii equation. The leading second order term is a 
perturbative Bogoliubov Hamiltonian 

£ P crt. = J dx (sft,-Sj>) C ( ) = 



E 



-n u n u n 



p2 

q 



+ 



2M q 2M 



(20) 



As is well known [6], this theory is predicting its own 
demise. The perturbative Hamiltonian (20) predicts in- 
definite spreading of phase and position with time, 

(9 2 ) ~ t,((q-q ) 2 ) ~ t, 

while at the same time the derivation of (20) requires 
them to remain small. In the next Section I introduce 
non-perturbative treatment of the zero modes that does 
not suffer from this inconsistency. 



V. NON-PERTURBATIVE THEORY 

The key point is observation that any classical field can 
be expanded as 
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Here 6„'s are complex Bogoliubov amplitudes. For the 
lower component to be a complex conjugate of the up- 
per component the coordinates jg and 7 g must be real, 
compare Eqs.(ll,15,16). The real collective coordinates 
9 and q are implicit in definitions of <f>o(x), adjoint modes 
and Bogoliubov modes. Unlike in the perturbative treat- 
ment (19), here 9 and q are not assumed to be small, but 
they can take non-perturbatively large values. 



A. Effective Hamiltonian 

The Gross-Pitaevskii equation (1) follows from a La- 
grangian 



l = J dx ( rt^M - |h^| 2 - gW + mH 2 



Substitution of Eq.(21) to L, expansion to second order 
in 7 and 6, and subsequent integration over x give an 
effective Lagrangian for the collective coordinates: 



hN9 + Pq - 



7 2 
iq 



2M 2M q 

Here N is a total number of atoms (conjugate to 9) 



(22) 



N=\i U\d <t>) - c.c] = N + ^ + 0(77, bj, bb) , 



(23) 

and P is center of mass momentum (conjugate to q) 
P= l -ih[{cl>\d q( j>)- c.c] = 7 , + 0(ry,67,W). (24) 



Legendre transformation leads to an effective Hamitonian 
H eS = Y e n b* n b n + JL + JL. (25) 



2M e 2M„ 

To complete the transformation 7's must be expressed 
as functions of the canonical momenta N and P: jg s=s 
h 2 (N — N a ) 2 and 7^ w P 2 . These approximate 79 and 
7 g result in a non-perturbative Bogoliubov Hamiltonian 

h 2 (N-No) 2 P 2 



2Mg 



2M„ 



(26) 



B. Quantum dark soliton 

H e s can be quantized by replacing c-numbers with op- 
erators. Non-zero commutators are: [6fe,6j] = S(k — p), 
[9, N] = i, and [q, P] = ih. With constant No we also 
have [9, N — No] = i. A quantum Hamiltonian in coordi- 



nate representation where N~N a 
is 

h 2 

H cS = £ tnbi b n - —dl 



-idg and P = —ihd„ 



dl . 



2M q q 



(27) 



This Hamiltonian is a sum of phonon terms, a phase dif- 
fusion term, and a soliton diffusion term with negative 
mass M q < 0. 
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In the framework of the non-perturbative theory of 
zero modes one is free to work in a subspace of Hilbert 
space with definite total number of atoms N. A state 
with definite N = N has a wave function ~ exp(iN9) 
which covers the whole range of 9 e [— n, +n) including 
non-perturbatively large values of 9. Such a wave func- 
tion, and consequently also a state with definite N, is 
beyond reach for the perturbative theory. In the sub- 
space with N — N the Hamiltonian becomes 



fjN 



2M q 9 



(28) 



In the same way, a wave packet for soliton position q can 
disperse covering most of the box without contradicting 
assumptions of the theory. For example, initial width 
Ag of a gaussian wave packet grows like 



Aq(t) = 4/Agg + 



h t 2 large* 



Ct 



4M„ 2 Aq 2 8 p Aq 



(29) 



This dispersion of soliton position becomes comparable 
to the soliton width, Aq = £, at the time Tdiffusion = 
8hAqo/g which depends on the initial dispersion Ago- In 
between t = and t — Tdiffusion a hole in single particle 
density distribution fills up with atoms, compare Fig. 2. 




\ % "R 

FIG. 2. Quantum diffusion of soliton position. This fig- 
ure shows the single particle density (ft (x)^(x)} / p at t = 
(dashed line) and at t = Tdiffusion (solid line) when the disper- 
sion of soliton position Aq equals the soliton width £. 



It turns out that there is non-zero minimal uncertainty 
of soliton position Ag m i n > 0. Suppose that we have 
two condensates with a dark soliton, but with different 
soliton positions gi and g 2 . What is the minimal distance 
between the solitons |gi — qi | when it becomes possible to 
distinguish these two condensates by a suitable quantum 
measurement? To answer this question we must calculate 
overlap between these two condensates of N atoms and 
see how it decays with the intersoliton distance, 



j>r qi \^r 92 ) 
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N 
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(30) 



Two condesates become orthogonal and in principle 
distinguishable when the intersoliton distance becomes 

greater than Ag m i n = • This is fundamental limita- 
tion derived only from properties of the quantum states 
and not of any particular measurement technique. This 
fundamental limitation means that the initial gaussian 
wave packet cannot be localized better than Ago = 
Ag m i n . This minimal dispersion leads to the minimal 
soliton diffusion time 



^diffusion — \ ~ 
C V 2 



(31) 



It is interesting to evaluate T m i n for the parameters of 
the Hannover experiment [2]. This condensate can be 
approximated by a quasi-lD harmonic trap like in Ref. 
[7] . In the present paper we consider uniform condensate 
in a box with linear density p. Using the linear density 
in the center of the effective quasi-lD trap as p we obtain 
Tmin = 8.0ms. This is close to the time O(10ms) when 
solitons are observed to loose contrast in that experiment. 

Finally, I make a brief comment on solitons mov- 
ing with finite velocity v with respect to condensate. 
Condensate wave function for x € {xl,xr} is 0o = 
e~' 9 VP (i@ + at'<mha x - q - vt ^ with = v/c and a = 

\J\ — 2 . Calculation of zero modes and adjoint modes 
follows the same lines as for v = 0. Mass of a moving soli- 
ton is Mg = aM q and the overlap between condensates 
decays when the intersoliton distance becomes compara- 
ble to Ag min = Ag m i n /a 3 / 2 . The minimal soliton diffu- 
sion time is 



^diffusion 



Tdiffusion 



(32) 



For Hannover solitons [2] moving with velocities v/c — 
0.4, . . . , 0.8 the minimal time is T min = 8.7 . . . 13.3ms. 



VI. CONCLUSION 

This paper develops non-perturbative theory of zero 
modes in the Bogoliubov theory of atomic BEC. In the 
non-perturbative approach phase fluctuations and fluctu- 
ations of soliton position are not restricted to be small. 
This theory predics that an initially well localized soli- 
ton wave packet is going to disperse beyond the soliton 
width. 
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When parameters of the present model are fit to the pa- 
rameters of Hannover experiment [2] , then dispersion of 
soliton position becomes comparable to the soliton width 
after 10ms from soliton creation. This number is consis- 
tent with earlier perturbative studies of soliton diffusion 
in a harmonic trap [7]. This diffusion time is also con- 
sistent with the time when the dark soliton appears to 
loose contrast or gray in Hannover experiment [2]. 

Interaction of the soliton with a thermal cloud was of- 
fered as an explanation [8] of the observed graying. The 
quantum diffusion described here and in Refs. [7] is a 
mechanism that operates even in the T = quantum 
limit when the thermal cloud is turned off. The quan- 
tum diffusion mechanism could be tested in an exper- 
iment with variable temperature. Another dissipation 
mechanism due to non-uniform condensate density was 
suggested in Rcf. [9]. The influence of the inhomogenity 
could be elliminated in an experiment with a static soli- 
ton (v = 0) in the center of the trap. In principle it is 
possible to make an experiment where the other mecha- 
nisms [8,9] are turned off leaving the quantum diffusion 
as the only known instability of the dark soliton. 
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